Thermal noise in superconducting quantum point-contacts 
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0\ . Abstract 



We present a theory for the frequency-dependent current fluctuations in 



^ ' superconducting quantum point-contacts (SQPC) within the dc transport 

o ■ 

^ ' regime. This theory is vahd for any barrier transparency between the tun- 

^ \ nel and balhstic hmits, yielding an analytical expression for the fluctuations 

I spectrum in the subgap region. It is shown that the level of noise in a quasi- 

ballistic SQPC may have a huge increase in comparison with the case of a 
normal contact carrying the same average current. The effect of this high 
level of noise on the actual observability of the current-phase relation for 
a ballistic point-contact is discussed in connection with recent experimental 
measurements. 
PACS numbers: 74.50.+r, 85.25.Cp, 73.20.Dx 

Present technologies make it possible to fabricate superconducting point-contacts in the 
nanometer scale. Examples of these kind of systems are the recently developed atomic size 
break-junctions and the split-gate superconductor - two-dimensional electron gas - su- 
perconductor junction of Takayanagi et al. 0. In both cases the electronic transport takes 



place through a reduced number of quantum channels, the contact transmission being a con- 
trollable quantity. These features make these systems very attractive for testing theoretical 
models of the superconducting transport beyond tunnel conditions. 

Recently, there have been a number of theoretical works devoted to a detailed analysis 
of both the dc and ac response of a single channel point-contact l^-^l. In particular, new 
illuminating results have been obtained for the ac current in the hitherto less understood 
limit of small bias voltages PJ^. However, little attention has been paid to the effect of 
thermal fluctuations in the transport properties of this kind of devices. This is an important 
issue both due to its intrinsic interest as a nonequilibrium phenomena and also because they 
limit the observability of the measured characteristics. It is clear that thermal fluctuations 
have to be taken into account if a direct comparison between theory and experimental 
results is to be carried out. Regarding this last point, some recent experiments have 
shown the deviation of the measured current-phase relation in a mechanically controllable 
break junction with respect to the theoretical predicted one 0. Some authors have recently 
pointed out the importance of thermal fluctuations as a source for this deviation 

The aim of this communication is to present a theory for the thermal current fluctua- 
tions of a SQPC in the dc regime valid for any contact transmission. This theory yields an 
analytical expression for the zero-frequency noise which, in the limit of low barrier trans- 
parencies, differs strongly from the standard tunnel theory result 0. In the opposite limit, 
i.e. for a ballistic contact, we find that the current fluctuations diverge when the supercur- 
rent tends to its maximum value. We claim that this fact explains the difficulties found 
for the experimental observation of the predicted current-phase relationship for a ballistic 
contact. 

In recent works we have introduced a theoretical approach for the study of the transport 



properties of superconducting nanoscale constrictions [^|10|- Iii this approach the system 
is described by a Hamiltonian written in a site representation, from which the microscopic 
Bogoliubov de Gennes equations can be derived |]10| . Within this model the normal transmis- 



sion coefficient through the constriction can be expressed in terms of microscopic parameters. 



allowing to establish a complete correspondence with other approaches based on scattering 
theory. 

For our present purpose of describing an atomic size contact, it will be sufficient to 
analyze the following Hamiltonian 

H = Hl + Hr + J2ite"^^'cljn^ + te-'^'^c^^CL.), (1) 

where Hl and Hr are the BCS Hamiltonians for the uncoupled electrodes (defined as L and 
R) , t is the hopping parameter which defines the normal transmission through the single 
quantum channel connecting both electrodes; and is the total superconducting phase 
difference between the electrodes. In the present calculations we shall neglect fluctuations 
in this superconducting phase difference and concentrate in the contribution to the current 
fluctuations arising from thermal excitation of quasiparticles. 

Within this model, the operator associated with the current through the contact can be 
written as 

where the different creation and annihilation operators appearing in Eq. (2) are the usual 
Heisenberg operators at a given time r. Then, the spectral density of the current fluctuations 
is defined as 

S{uj) = hj dre'"^^ [< (5/(r)5/(0) > + < (5/(0)5/(r) >] , (3) 

where (5/(r) = /(r)- < / >. 

For the evaluation of the above averages, we perform a decoupling procedure which is 
consistent with the BCS mean field theory. The spectrum S{u!) can then be expressed 



in terms of the single-particle non-equilibrium Green functions [^, G^'^ (fi) and G^j^ 
(where a and (3 can be either L or R). In a superconducting broken symmetry representation 
ig Gtpi^) is defined by 

G+-(^])= /rfre^^^G+-(r,0), (4) 



with 



G'+«-(r,0) 



[<cl^iO)ci^ir)> <c^i(0)4( 



r > 



and G-+(r,0) = [G+-(0,r)]t. 

In terms of the functions G'^j^{Q) and G~j^{Q), S{uj) adopts the form 

^2 



S{u) = jjdn Tr[tG^l{n)G],tin + u;)i + 

iGtj,in)G-^U^ + uj)i - iGt^in)Gj,u^ + uj)i - 
iG+],{n)Gitin + iu)i + idj^nn + u;)G],t{^)i + 
idtj^in + uj)G-LU^)i - iGl-^in + uj)G-^j,{n)i - 

iGtR{^ + u)G-^t{Q)i], 



(5) 



where t is the hopping interaction between the electrodes written in the (2 x 2) Nambu 
representation 
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te'^ 
-te^'^ 



(6) 



In the present communication we concentrate in the zero voltage case in which the average 
current is due to Cooper pairs. For the calculation of the nonequilibrium Green functions 



appearing in Eq. (5) we can the use the relations [10 



a, 13 



(Q) = Gun) - Gun) f{n) 



= - Gi^{n) - G:,p{Q) [i-fm 



(7) 



where f{^) is the Fermi factor and G'^j^'' are the retarded (advanced) Green functions of the 
coupled contact. These last quantities can be obtained up to infinite order in the coupling 
parameter t by solving the following Dyson equation 
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where S^i"^ — ^^rr — '^^lr — (^/jl"*) — ^- The indexes a, (3 and 7 can be either L 
or i?, and g'^^a^ are the retarded (advanced) Green functions corresponding to the left and 
right uncoupled electrodes. 

For the symmetric case, both electrodes have the same modulus of the superconducting 
order parameter, A, and these Green functions can be expressed as 

^LL M = 9rr\^) = 1 _ , = X 



\ 

(9) 



—UO ±17] A 



A —uj ±11] J 

where W is an energy scale related to the normal density of states at the Fermi level by 
p{^f) = 1/ {tcW) and 77 is a small energy relaxation rate that takes into account the damping 
of the quasi-particle states due to inelastic processes inside the electrodes. This parameter 
can be estimated from the electron-phonon interaction to be a small fraction of A ||13|. It 
is useful to define the normal transmission coefficient of the contact, which in terms of W 
and t has the form a = A{2t/Wf/{1 + {2t/WYY The spectral densities that are 



obtained from Eq. (8) are no longer singular at the gap edges and exhibit poles inside the 



superconducting gap, located at energies ujs = ±Awl — a sin (0/2), corresponding to the 



interface bound states |T5|. As stated in previous works, these bound states carry all the 
Josephson current in the limit of a short constriction [|10],|16|. Therefore, their contribution 
to the zero-voltage current fluctuations can be expected to be crucial, as is certainly found. 

Once the single particle Green functions are known, the spectrum S{uj) can be calculated 
using Eq. (5). The typical form of this spectrum is illustrated in Fig. 1, where S{uj) is 
plotted for fixed temperature and three different contact transmissions. Notice that for 
uj < 2A the spectrum is formed by two resonant peaks at a; = and uj = 2uJs, arising from 
the existence of the bound states at uJs- Qualitatively, the peak at zero frequency increases 
with increasing transmission, while the one at 2ujs is negligible for both nearly perfect and 
very small transmissions, adopting its maximum value around a ~ 2/3. For uj > A + \ujs\ 
contributions from the continuous part of the single particle spectrum become important. 



In the limit of a very weakly damped contact, i.e. t] ^ aA, it is possible to evaluate 
S{lj) at o; = and u = 2ljs analytically. We find 

^, , 2e^ vr A^a^ sin^(</)) . , , 

-5(0) = -T-- -^^/(^5) [1 - /M] , (10) 

At 1^ S 

and 

Qro ^ 2e2vrAVa-a)siiA|) r 
^(2-5) = ^- [/K) + fi-u^s) 

These expressions clearly display the important role played by the interface bound states 
in fixing the magnitude of the current fluctuations for subgap frequencies. It should be 
stressed that, although the absolute size of the current fluctuations depend on the estimated 
value of parameter r], its precise variation with the superconducting phase difference and 
temperature is controlled only by the contact transmission a. 

Our analytical results are strictly valid in the limit r] <^ aA and differ strongly from 
the equilibrium fluctuations obtained using standard tunnel theory [^], which yields S{0) ~ 
a(l + cos(0)) In A/r]. This last expression becomes accurate just in the opposite limit, 
rj ^ a A which holds in the tunnel regime, i.e. a -C 1. In ref. we have explicitly shown 
that the limits r] and a do not commute. This behavior can be understood in 
the following way: when rj <^ aA, multiple Andreev scattering processes give the dominant 
contribution to any dynamical quantity and should be included up to infinite order. On 
the other hand, when a is small enough (in such a way that aA <^ 77) these high order 
scattering events become heavily damped and the lowest term of the perturbative expansion 
in t gives the correct result. For a realistic SQPC in which, as commented above, 77 can be 
estimated to be a small fraction of A, the situation would always correspond to the weakly 
damped regime, except for extremely small values of a and therefore Eqs. (10) and (11) will 
accurately describe the low frequency noise. 

The analysis of Eq. (10) reveals some remarkable physical consequences. To begin with, 
and in contrast to the normal case where a reduction of noise is found |l^, S{0) experiences 
a dramatic increase when approaching the ballistic regime. More precisely, there is a value 



(11) 



of a, given roughly by the condition ksT ~ Av^l — a, above which there is an exponential 
increase of the thermal noise. On the other hand, in this last situation, there appears a very 
strong asymmetry on the phase- dependence of S'(O), with its maximum value progressively 
moving from = 7r/2to0 = 7r. This remarkable behavior should certainly have implications 
in the actual observability of the supercurrent-phase relation in a SQPC 

In order to analyze the importance of these thermal fluctuations it is convenient to study 
the ratio S{0)/2e < I{(f)) >, where < > is the phase- dependent average supercurrent 
(let us recall that for a normal contact the classical shot noise limit corresponds to 5(0) = 
2e < / >), given by |[lO| 



?/ ,N CTT A^a sin(0) , J cje L 

<'<^)>=T^;7r'"'"y^'- '''' 

In Fig. 2 we plot S'(0)/2e < /(0) > as a function of the superconducting phase difference 
for increasing values of the transmission a and two different temperatures. This figure 
illustrates the huge increase of thermal noise when the transmission becomes sufficiently 
large. As can be observed, for a reasonable choice of parameter t] and depending on the 
temperature, the level of noise can reach values several orders of magnitude larger than 
2e < /(0) >. When lowering the temperature this level of noise is reduced, but it will 
always be significant close to the ballistic case in a phase interval around = tt, just in the 
zone where the average current has its maximum at low temperatures. In fact, taking the 
limit a 1, the ratio 5'(0)/2e < /(0) > has the form 

(„ ^ 1, . H^^l^lSl , (13) 



which clearly diverges when approaches tt. The fact that the zero frequency noise has 
large values in the zone where the maximum of the average current occurs can explain the 
experimental difficulties found to observe the predicted ~ sin 0/2 form of the current-phase 
relation for junctions with direct conductivity, as reported in ref. [^. 

Finally, it is worth discussing our result for the zero frequency noise of a weakly damped 
contact in the light of the Callen-Welton fluctuation- dissipation theorem. In general, this 



theorem relates the equihbrium current fluctuations with the linear conductance G by S{0) = 
ikBTG. This relation allows us to calculate in a straightforward way the phase- dependent 
linear conductance of a SQPC from Eq. (10) 



h ksTr] 



A2asin(0) us 
■sech(- 



(14) 



Eq. (14) coincides exactly with the result of ref. P] in which a direct calculation of the 
linear conductance of a SQPC was performed (this expression for G(0) has been recently 



rederived in [|l^] for the particular case of a ballistic contact). 

In conclusion, we have developed a theory of the thermal fluctuations for a SQPC in the 
dc regime. The noise spectrum exhibits resonant peaks at subgap frequencies associated with 
the existence of bound states in the constriction region. For the case of a weakly damped 
contact {rj <^ aA), we have obtained a closed analytical expression for the weight of these 
resonant peaks. We have shown that a striking consequence of the presence of these bound 
states is a huge increase of the low frequency noise level when approaching the ballistic limit, 
this high level of noise being particularly important when the average supercurrent is close 
to its maximum value. We claim that these results may explain the reported difficulties in 
measuring the predicted sin 0/2 behavior for < /(0) > in the case of highly transmissive 
contacts. Finally, we have discussed the connection between the present theory and the 
phase-dependent conductance of a SQPC by means of the fluctuation-dissipation theorem. 
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FIGURES 

FIG. 1. Current fluctuation spectrum of a SQPC in the dc regime for three different values of 
the transmission. The superconducting phase difference corresponds in each case to the maximum 
supercurrent and the temperature is ksT = 0.2 A. 

FIG. 2. The ratio between the zero frequency noise and the average supercurrent as a function 
of the superconducting phase difference for increasing values of the transmission coefficient. 
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